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Abstract. In this talk we will present the solution of the two-channel Anderson impurity 
model, proposed in the context of the heavy fermion compound UBcia, and discuss 
briefly the more general multi-channel case. We will show results for the thermodynamics 
in the full range of temperature and fields and make the connection with the current 
experimental situation. 



1. Models of Uranium ions 

During the last two decades, there has been growing interest in materials 
whose desciption falls outside the framework of Landau's Fermi Liquid 
Theory. Examples are the high Tc superconductors, heavy fermions and 
quasi- ID conductors. 

Uranium heavy fermion compounds provide very interesting examples, 
exhibiting unusual specific heat temperature dependence, 7 = Cy/T ~ 
InT. This behavior is often related to multichannel Kondo physics, occa- 
sioned in the case of Uranium by the ground state configuration which is 
degenerate and non-magnetic, so that one is lead to consider a quadrupolar 
Kondo scenario (1, 2). 

The starting point is to consider the energies of the uranium ion in 
its different states of valence. They fall on a parabola depicted in Fig. 1. 
The ground state corresponds to a U^"*" ionization state, with a Sf'^ shell- 
configuration. Considering the spin-orbit Hund's coupling, the lowest mul- 
tiplet corresponds to J = 4. Taking, further, into account the splitting due 
to a cubic crystal field (as is the case for the heavy fermion compound 
UBeis) one finds that the lowest multiplet corresponds to a non-Kramers 
Fa doublet (see Fig. 1). 

A similar analysis for the neighboring states of valence indicates that 
they corresponds to: (i) U^+: a 5f^ shell-configuration that splits giving a 
J = 9/2 lowest multiplet configuration, that further splits giving a lowest Fg 
Kramer's doublet and (ii) U^^: a 5f^ shell-configuration that splits giving 
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Figure 1. Schematics of the energy and lowest angular momentum multiplet for different 
states of valence of a uranium ion. 




Figure 2. Schematic depiction of the spin-orbit and crystal field splittings for the ground 
state valence configuration of uranium in a cubic medium. 



a J = 5/2 lowest multiplet configuration, that further splits giving a lowest 
T^ Kramer's doublet. These three valence states are in fact very close in 
energy and their relative positions are not completely resolved, since the 
experimental evidence is not conclusive and sometimes contradictory (3, 4). 
The most accepted scenario is that of a T-^ quadrupolar doublet ground state 
coming from the tetravalent state and next a Fg magnetic doublet coming 
from the trivalent state of uranium. Neglecting the contributions from the 
pentavalent state of the uranium as well as those from excited crystalline 
electric field states we can write down the following Hamiltonian, the two- 
channel Anderson impurity model (see Fig. 1 for a schematic depiction). 



a = +. 



Figure 3. Schematics of the electronic configurations and symmetry multiplets retained 
for the Hamiltonian of UBeia. 
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With Q = 5]6Ua + E/^/- = l 

a cr 

where the bar on top of the index a indicates that it transforms according 
to the complex conjugate representation. We used slave-boson language and 
denoted Sq = Er^ and £5 = (let us also define e = Sg — Sq and A = 
V'^/2). Notice that the model has full SU{2) SU{2) internal symmetry. 
In the parameter regimes when |e| S> A, we can map it via a Schrieffer- 
Wolff transformation into a quadrupolar two-channel Kondo model (with 
a localized quadrupolar moment in the case when Sq <C £s), or a magnetic 
two-channel Kondo model (with a localized magnetic moment in the case 
when Eq 3> £s), in both cases with a coupling constant given by J = 
/ \£s — £q\ (assuming zero chemical potential). 

Wc will also consider a generalized version of the Hamiltonian, the multi- 
channel Anderson impurity model. In slave-boson language the Hamiltonian 
looks exactly as in the two-channel case, but this time the indices take 
values according to: (i) cr = 1, A?^ and (ii) a = 1, . . . , M. If we follow 
the standard nomenclature wc would call a the spin index and a the channel 
index, and the model will be the SU (N) Cg) SU (M) version of the multi- 
channel Anderson model (5). 

As in the two-channel case we can perform Schrieffer- Wolff transforma- 
tions to derive effective models in any of the two local moment regimes 
of the problem. When <^ Eg (assuming always zero electronic chemi- 
cal potential) we obtain an 'iV-channel SU{M) Coqblin-Schrieffer model', 
whereas when Eg ^ Eq we obtain an 'M-channel SU{N) Coqblin-Schrieffer 
model'. 

This model has been extensively studied over the past fifteen years by 
a variety of methods: Large-N methods (6), NCA (4), conserving T-matrix 
approximation (7), Monte Carlo (4) and Numerical Renormalization Group 
(8) among others. All these techniques have difficulties accessing the mixed 
valence regime. 

Wc showed that this model is integrable and carried out a full analysis 
by means of the Bethe Ansatz technique. That allowed us to obtain exact 
results characterizing the ground state of the model as a function of its pa- 
rameters and very accurate numerical results for the thermodynamics of the 
model for all temperatures and fields, as well as all valence regimes, includ- 
ing the mixed valence. We stress that is in the mixed valence regime where 
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the Bethe Ansatz solution is most crucial, since all the other techniques 
have difficulties accessing it. 

2. The Thermodynamics of the Model 

In this section we present the results of the numerical solution of the 

Thermodynamic Bcthc Ansatz equations for the two-channel case (i.e. N = 
M = 2). This is the first model proposed in the context of UBci3(l). We will 
study the physics of the model on its own right and also briefiy discuss its 
virtues and inadequacies to describe the physics of the uranium compound 
that motivated it. The conclusion that we will reach is that one needs to 
go beyond the simplest scenario given by the two-channel model. 

In Fig. 2 we show the behavior of the impurity entropy as a function of 
temperature for different values of (e — /x) , and in Fig. 2 the reader can see 
the effect of switching on an external field. 
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Figure 4- Impurity contribution to the entropy at zero field and as a function of 
temperature for different values of e — /i. As the temperature goes to zero all curves 
approach the universal value ks In v^. Positive and negative values of e — /^ fall on top 
of each other. 

At high temperatures the entropy is 5imp = fee In 4 in agreement with 
the size of the impurity Hilbert space. For |£ — /u| 2> A the impurity entropy 
is quenched in two stages. The degrees of freedom corresponding to the 
higher energy doublet are frozen first. The entropy becomes S'imp = fee In 2 
and the system is in a localized magnetic or quadrupolar moment regime 
depending on the sign of (e — ji). As the temperature is further decreased, 
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Figure 5. Impurity contribution to the entropy in the presence of an applied quadrupolar 
field and as a function of temperature for different values of e — As the temperature 
goes to zero all curves approach zero. Positive and negative values of e — /i fall on top of 
each other only at high temperatures and, asymptotically, for low temperatures. 

the remaining degrees of freedom undergo frustrated screening leading to 
entropy Simp = A;Bln-\/2. On the other hand, for values of le — /x| <^ A, 
the quenching process takes place in a single stage. The initial stage of the 
quenching process happens when the system makes the transition from a 
high-energy state of valence = ^ (i.e. when both states of valence are 
equally likely) to the state of valence that will be present at low-energies, 
n^(e — /Lt). This is further illustrated in Fig. 2. 

As (e — ji) is varied, the behavior interpolates continuously between 
the magnetic (n° = 1) and the quadrupolar (n° = 0) scenarios. The zero 
temperature entropy is found to be independent of e in accordance with 
our analytic results (9). External fields, either magnetic or quadrupolar, 
constitute relevant perturbations that drive the system to a Fermi Liquid 
fixed point with zero entropy. 

In Fig. 2 we show the impurity contribution to the specific heat. The 
two stage quenching process gives rise to two distinct peaks. The lower 
temperature peak is the Kondo contribution centered around a tempera- 
ture Ti {e — n), whereas the higher temperature peak - often referred to 
as the Schottky anomaly - is centered around a temperature (e — /x). 
Approximate expressions for T/j ^ can be read off from the curves: 
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Figure 6. Charge content of the impurity site as a function of temperature for different 
values oi e — fj,. At high temperatures all curves start from ric = 1/2 and flow, as the 
temperature is lowered to speciflc values n° {e — ji). 



with 1 < a < 4. For large |£ — /i| the two peaks are clearly separated and 
the area} under the Kondo peak is ks ln\/2 while that under the Schottky 
peak is In 2. 

As mentioned earlier, the model was proposed as a description for 
the uranium ion physics of UBeia. It is expected to describe the lattice 
above some coherence temperature. We provide in the inset of Fig. 2 the 
experimental data for the 5f-derivcd specific heat of the compound. It is 
obtained by subtracting from its total specific heat, the specific heat of the 
isostructural compound ThBci3 containing no 5f electrons (10, 11). This 
way one is subtracting the phonon contribution as well as the electronic 
contribution from electrons in s, p and d shells (the procedure is quite 
involved and we refer the reader to the cited articles for full details). The 
sharp feature at ~ O.Sif signals the superconducting transition of UBeia 
and falls outside the range where this compound might be described by 
considering a single impurity model. 

In order to be certain of having eliminated lattice effects it would be 
better to carry out the measurements on Ui-xThj^Beis. For x > 0.1 the 
compound has no longer a superconducting transition and the lattice coher- 
ence effects are largely suppressed. Further, there are several experimental 



^ Here by area we mean the following integral: J Cv {T) d(lnT). In other words, the 
area as plotted, i.e. in a logarithmic scale. 



7 




Figure 7. Main plot: impurity contribution to the specific heat as a function of tem- 
perature for different values of e — As this parameter approaches zero the Kondo 
contribution (left) and the Schottky anomaly (right) collapse in a single peak. Inset: 
experimental data for the 5/-derived specific heat of UBeis . 



indications that support the idea of an impurity model description of the 
thoriated compound for a wide range of temperatures (12). 

Concentrating on the temperature range containing the Kondo and 
Schottky peaks wc conclude that no values of £ and A of the 2-channel 
model yield a good fit. Further, the entropy obtained by integrating the 
weight under the experimental curve falls between fee In 4 and /cb In 6. This 
suggests that a full description of the impurity may involve another high 
energy multiplet (possibly a r4 triplet, cf. with the work of Koga and Cox 
(8)) close to the Fs to yield the peak for the Schottky anomaly, with the 
Tq doublet falling outside the range of measurements. The nature of the 
multiplet could be deduced from further specific heat measurements. For 
an n-plet degenerate with the F3, one has an SU{2) SU{n + 2) Anderson 
model and the area under Cl^^/T is then given by ln[^^^ sec ^^], while 
if the n-plet is slightly split off the doublet, the area is given by fee ln[^^] 

(13). Also the scale of energy involved is then determined by crystal field 
splitting and corresponds to the scale observed. The split off non-magnetic 
triplet (i.e. n = 3 in the configuration seems, indeed, to be the most 
reasonable scenario for UBei^. We are currently working in this direction. 
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